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Minimal conductivity in graphene: interaction corrections and ul- 
traviolet anomaly 
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Abstract. - Conductivity of a disorder-free intrinsic graphene is studied to the first order in the 
long-range Coulomb interaction and is found to be a — <to(1 + 0.01<?), where g is the dimensionless 
("fine structure") coupling constant. The calculations are performed using three different methods: 
i) electron polarization function, ii) Kubo formula for the conductivity, iii) quantum transport 
equation. Surprisingly, these methods yield different results unless a proper ultraviolet cut-off 
procedure is implemented, which requires that the interaction potential in the effective Dirac 
Hamiltonian is cut-off at small distances (large momenta). 



Introduction. — Low-frequency optical conductivity 
of undoped (intrinsic) graphene free of disorder is known 
to have a universal value of cto = e 2 /4h~ [1-13]. Experi- 
mental measurements [14,15], which yielded a value some- 
what bigger than the theoretical predictions, motivated 
the studies of the possible role played by electron-electron 
interactions. The findings of Ref. [16] that the combined 
effect of self energy (velocity renormalization) and vertex 
corrections leads to a suppression of the optical conductiv- 
ity at low frequencies have been questioned in Refs. [17, 18] 
on the basis of scaling arguments. The latter indicate 
that the large logarithmic (momentum cut-off dependent) 
terms in the self-energy and vertex corrections cancel each 
other. We note that Ref. [16] and Refs. [17, 18] agree on 
this cancellation in the lowest order in electron-electron 
interaction but differ on whether the higher order terms 
feature similar cancellation. It appears that the analysis 
of Ref. [16], though valid in the first order, fails for higher 
orders, and that the conclusion of the suppression of the 
conductivity at low frequencies is not valid. 

The theory presented in Refs. [17, 18] implies that the 
low-frequency dependence is properly described by the 
lowest order correction. Indeed, to the first order in in- 
teraction the conductivity is expected to yield, cr/<To = 
1+Cg, where C is some constant, g — e 2 / nv is the interac- 
tion strength; k is the dielectric constant of a substrate and 
v is the electron velocity in graphene. Renormalization 
group approach for 2D Dirac fermions predicts that the 
interaction strength g is a running coupling constant that 



depends on frequency g — ► g(u) [20,21]. At low frequen- 
cies g((j) flows to zero, so that higher order corrections to 
the electron velocity become progressively negligible and 
it is sufficient to consider only the first order renormal- 
ization of velocity (electric charge is not renormalized) : 
g(u>) = g/[l + | In QCv/u)], where K is the momentum 
cut-off. Combining these expressions gives, 



a/a = 1 + 



Cg 



1 



In {Kv/ui) ' 



(1) 



with the low-frequency behavior of the conductivity being 
determined by the constant C alone. Calculation of this 
constant, therefore, becomes an important task. While 
Ref. [17] did not calculate C, Ref. [18] provided the fol- 
lowing value 

C= Hz*E„ .51. (2) 
12 v ; 

This result predicts quite a considerable variation of a 
with the frequency for typical values of the bare graphene 
interaction constant g (which can exceed 1). 

In the present Letter we test the above prediction (2) by 
performing a perturbative calculation of the minimal con- 
ductivity to the first order in electron-electron interaction 
using three different methods, based on, a) electron po- 
larization operator, b) Kubo formula for the conductivity, 
c) kinetic equation. We point out that crucial anomaly, 
which does not appear in a non-interacting case, occurs for 
the interaction correction. Three above mentioned meth- 
ods would give essentially different values for the constant 
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Fig. 1: Self-energy, a), and vertex correction, b), to the con- 
ductivity <t(w) and polarization operator Yl(uj,q) in the first 
order in electron-electron interaction (dashed line) . The vertex 
(black dot) is equal to 1 in the case of the polarization operator 
and to ev& in case of the conductivity. The two quantities are 
related to each other by the particle conservation condition, 
Eq. (5). 



C unless some appropriate high-momentum cut-off proce- 
dure is implemented. We argue that expression (2) over- 
estimates the interaction correction by almost two orders 
of magnitude and show that the numerical value of C is 



C = 



19- 6tt 
12 



0.01 



(3) 



We will now proceed to demonstrate that the difference 
between Eqs. (2) and (3) originate from handling of sin- 
gular integrals at large electron momenta. 

The first method to be presented is based on the calcu- 
lation of electron polarization operator and has an advan- 
tage of being free from any such singular integrals. 

Polarization operator. — Single intrinsic 2D 
graphene layer is described by the chiral Hamiltonian, 

fl = «E^-p^E ^ + <M£>> (4) 

ip ijpkq 

where "hats" denote operators in a pseudo-spin space (<x 
represents the usual set of Pauli matrices), the sum over 
Latin indices is taken over two nodal points and two (true) 
spin directions. The interaction potential is Vq = 2ne 2 j 'nq; 
we also denote, ^ p = J cPp/(2n) 2 7 and set K=l. 

First-order interaction corrections to the conductivity 
are given by the two diagrams shown in Fig. 1, with the 
vertices denoting the operators of electric current, ev&. 
Another possible method to derive the homogeneous opti- 
cal conductivity is to calculate the corresponding diagrams 
for the electron polarization operator H(u>, q) and then uti- 
lize the particle conservation condition, 



a(u) = lim ^-^-U(w,q). 
g^a q 2 



(5) 



The calculation of the polarization operator to the first 
order in g requires two diagrams [19], 



pp' 



2G e pG e 'p'G ep 



XGe-l-^.p+q -+- G £ pG e 'p'G e '-f w .p'-|-qG c + Wj p+q 



(6) 



where the electron Green's function in the subband repre- 
sentation is 



G™ = - V ' ~" >ap 



2 e - (3(vp - in) ' 



(7) 



with <7 p = & ■ rip being the projection of the pseudo-spin 
operator onto the direction of the electron momentum 
n p = p/p. Factor 4 in Eq. (6) accounts for the (true) 
spin and valley degeneracy. Taking energy integrals and 
performing pseudospin trace operation we obtain for the 
first term in Eq. (6), 



U a (uj, q) = 2 ^ ^p-p' cos0 PiP / 



1 — cos 6 



p,p+q 



pp'P 



{uj + 2(3vp - i(3n) 2 



■ (8) 



In the expression (8) we kept only terms that lead to 
the lowest order contribution in the external momentum, 
Il a (u),q) cx q 2 , which are necessary for the calculation of 
the homogeneous conductivity. Using Eq. (5) we obtain 
the corresponding contribution, 

2_, ^p-p'COS0 p , p y^ 2 _ 4 ^ 2)2 , (9) 

where the frequency is presumed to have a positive in- 
finitesimal imaginary part. The second term in Eq. (6) is 
evaluated similarly, 



a a (uj) = ie 2 u J ' 



U b (u,q) = -2J2 



m 2 (n p -n p+q ) ■ (iy - iy +q ) 

(oj 2 - 4v 2 p 2 )(ui 2 - 4v 2 p' 2 ) 



Av 2 pp' sin flp, p+q sin fl p ', p ' +q 

(lu 2 - Av 2 p 2 )(uj 2 - Av 2 p' 2 ) 



W- (io) 



Expanding to the quadratic order in q we obtain the vertex 
correction, 

2-t cos 2 9 P p > + Av 2 cos P p > 

To the first order in interaction the conductivity is given 
by the sum a = a + a a + a b . The second term here, given 
by Eq. (9), contains a strong divergence at p = u/2v. This 
divergence, however, is simply a consequence of the renor- 
malization of the electron velocity by electron-electron in- 
teractions. To make the integrals regular we note that 
both the zeroth-order term [19] and a a can be written as 



(To + a a 2ie u) ^ ^ ^ 2 _ 4 v 2p2) ' 



(12) 



with v p = v + ± J2 P > V p-p' cos ^p,p' = v[l + f In (K,/p)} 
being the renormalized velocity (where K. is the upper mo- 
mentum cut-off). Indeed, expanding the integrand to the 
first order in the interaction one recovers Eq. (9). Note 
that the value of v p coincides with the electron velocity 
found from the perturbation expansion for the electron 
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Green's function [20]. The integral in Eq. (12) is regular. 
Calculating the real part of Eq. (12) we obtain 



°0 + CT„ = (T (l + 



(13) 



Note that the interaction correction in Eq. (13) is due to 
the curvature of electron spectrum. 

Calculation of the real part of Eq. (11) can be reduced 
to the following dimcnsionless integral (x — 2vp/u), 



Ob = -°o .9 



f dO f da 
J ~ J (1-x 



dx cos 9(x + cos ( 



*)y/x 2 + 1 - 2x cos9 



(14) 



where the integral is taken in the principal value sense. 
Using the identity 



dx(x + cos ( 



_ sin 2 (fl/2) 

(1 - x 2 )Vx 2 + 1 - 2xcos(9 ~~ cos (6/2) 



and integrating over the angle 9 we obtain, 

8 - 3tt 



°b = (T g- 



ln [tan(0/4)], 
(15) 



(16) 



Combining Eqs. (13) and (16) we finally arrive at Eq. (3). 
The vertex correction is negative and nearly cancels the 
self energy correction. The frequency independence of 
Eqs. (13) and (16) is analogous to the independence of 
the non-interacting conductivity <jo- 

Kubo formula for conductivity. An advantage of 
deriving minimal conductivity from the polarization oper- 
ator originates from the fact that large logarithmic contri- 
butions do not appear in different terms in this formalism. 
On the other hand one could begin with a straightforward 
application of the Kubo formula, which as we will see, 
does not offer such a simplification. As a result one has 
to deal with logarithmic contributions which ultimately 
cancel. The starting expression is the expression for the 
optical conductivity 



K(u) - K(0) 



(17) 



via the current-current correlation function, which is given 
in the zeroth order by 

/j 
^- a x G ep a x G e+UiP , (18) 

and in the first order by 

K a {u) = 8*eVTr W ^| W 
pp' 

XO'jGepGe'p'O'jGtpGe+u.p, (19) 

K b { U ) = 4*eVTrW|^ W 
pp' 

x o x G €p G e i p ' a x G £ i +UJ p ' G e+LU p , (20) 



here the subscripts a, b denote the contributions from the 
self energy and vertex diagrams, Fig. 1. Note, however, 
that the corresponding contributions into the conductiv- 
ity, which we denote here by a a and a b , do not satisfy 
the condition (5) term by term. However their sum has to 
obey it, <r a + a b = a a + a b . 

Calculation of Eqs. (18-20) is similar to the above 
derivation for the polarization operator. 



(T a (w) 



ie 2 ui Vp-p' cos t 



12v 2 p 



2^2 



U) 



PP' 



P 2 (uj 2 — Av 2 p 2 ) 2 ' 



(21) 



a b (u>) 



ie 2 ui 



p' 



cosf 



-^(lu 2 -8v 2 p 2 )cos6 -4v 2 



(ui 2 - Av 2 p 2 )(ijj 2 - Av 2 p' 2 ) ' 

I' I' 

(22) 

where we omit the subscripts in cos0 P! p/. The expres- 
sions (21,22) are to be contrasted with Eqs. (9) and (11). 
The obvious distinction arises from the fact that the in- 
tegrals in Eqs. (21) and (22) are logarithmically divergent 
though these divergencies ultimately cancel in their sum, 
<r(w) = a a (uj) +(Tb(w)- However, a more striking obser- 
vation can be made if one calculates the difference of the 
two expressions, 



(t(uj) — a(ui) = I, 



where 



p — p cos t 



I = 2ie us y K,-p' cos 0—77—0 — 



PP 



p 2 p'(uj 2 — 4v 2 p 2 ) 



(23) 



(24) 



Before addressing the issue of a numerical value of / let 
us briefly describe the third method for the calculation of 
the conductivity. 

Kinetic equation. — The kinetic equation to the low- 
est order in electron-electron interaction and in the pres- 
ence of electric filed has the form [16] 



+ ivp[*P, /p] + eE ' |^ ^ 1 E V p-p' t/p' > /p1' ( 25 ) 

p' 



dp 



where / p is the 2x2 matrix distribution function. The sec- 
ond term in the left-hand side represents the rate of change 
of the electron distribution function during its preces- 
sion in the momentum-dependent "pscudo-Zeeman" field. 
The third term is the usual drift in the momentum space 
caused by external electric field. Finally, the right-hand 
side account for the exchange electron-electron interac- 
tion (Hartree contribution being zero by virtue of electric 
neutrality and spatial homogeneity). Given the solution 
of kinetic equation (25) one can find electric current and 
optical conductivity from j = ct(lj)F, = 4eiTr^ p / p <r. 
The detailed solution of Eq. (25) to the first order in E 
and Vp-p' was found in Ref. [16]. For the conductivity it 
yields, 



(Thin M = Sie 2 u; ^2 v p-p' cos # 7 2 \ 



(lu 2 - Av 2 p 2 ) 2 
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_ A - 2 v 2 pcos 2 6 + v 2 p' cos 6 

This expression is different from both <t(uj) and a(ui). In- 
terestingly, 

a{u)-a kin {u) = I/2. (27) 

Discussion. — Three different values obtained from 
the polarization operator, o~(u>), Kubo formula, u{uS), and 
kinetic equation, <Thi n {u), respectively, indicate an incon- 
sistency of the theory of interacting two-dimensional Dirac 
fcrmions unless 1 = 0. We will now demonstrate that the 
conclusion of whether I = or I ^ depends on the way 
the ultraviolet cut-off is imposed in the calculation of a 
singular integral over p' in Eq. (24). 

(i) Hard cut-off. Let us first assume that the divergent 
momentum integral is extended only to p' < JC. By noting 
that interaction potential depends only on s = (p — p') 2 
and that p' — pcos9 = jds/dp' we then obtain for the 
latter integral in case when V(s) — 2ire 2 /^/s, 

2ir K 

Em Q P~PCos6 e 2 f de f dp' ds 

Vp-p'COS# = — / — COS& / — — — — - 

, p p p' 2 J 2tt J yfidp' 

P 

= e 2 / ^ cos V^K^pj. (28) 
J 2tt 
o 

Expanding \J s(JC,p) w JC — p cos 9 for large values of JC, we 
observe that the integral here is cut-off independent and 
equals —e 2 p/2. It is now straightforward to verify that 
equation (24) gives I = ~<T g/2. Such value of I yields 
Eq. (2) reproducing the result of Ref. [18], and precisely 
accounts for the difference between Eq. (2) and our result 
(3). However, as shown above, such an ultraviolet cut-off 
yields three different values of the conductivity depending 
on which method is being used and is therefore unphysical. 

(ii) Soft cut-off. The anomaly encountered in the ex- 
pression (24) is specific for V q cx q^ 1 behavior of the in- 
teraction potential. For any faster decay of interaction at 
large momenta the integral I vanishes. Let us demonstrate 
this point by assuming 

2-7TP 2 

V(q) = —e-i' K , JC -> oo. (29) 
Calculation similar to the preceding one gives, 

2-7T OO 

e 2 f d8 f e -vWC Q s 

lim - / ^cos0 / dp'f__|i = 0, (30) 
K^oo 2 J 2tt J Vs dp' K J 



so that 1 = and all three methods yield the same value 
(3). Similar conclusion will be reached if one assumes 
V(q) cx q~ 1 ~ v and subsequently takes the limit ?y — > 0. 

Having established that the hard cut-off utilized in 
Ref. [18] in the course of Kubo calculations actually re- 
sults in different and hence inconsistent results when other 



methods are used, it is time now to discuss the origin of 
this inconsistency. Terminating momentum integrals at 
some value p = JC means in fact an essential modification 
of electron spectrum at large momenta that effectively ex- 
cludes these states from possible virtual processes. Such a 
procedure, though not necessarily incorrect, can be made 
self-consistent only if it is accompanied by the appropri- 
ate change in the operators of electric current. Otherwise, 
the Ward identity, which ensures particle conservation, is 
violated. This is why the polarization function method, 
which does not involve current vertices in the course of 
calculations, gives results (in the form of convergent inte- 
grals) independent of the cut-off procedure. On the other 
hand both the Kubo formalism and kinetic equation do in- 
volve current operators explicitly and thus fail if the hard 
cut-off is implemented without a proper modification of 
current vertices. 

To the contrary, the soft cut-off procedure presented 
in this Letter does not require modifications of the elec- 
tron spectra (Green's functions) nor of the electric cur- 
rent vertices. It is thus self-consistent and quite naturally 
yields identical values for the conductivity irrespective of 
the method used. 

Conclusion. — We have calculated the first order 
interaction correction to the conductivity of intrinsic 
graphene. Within the Kubo and kinetic equation for- 
malisms the self-energy and vertex corrections contain 
large logarithmic frequency-dependent terms which ulti- 
mately cancel each other. Within the more convenient 
approach based on the calculation of the polarization op- 
erator, such terms do not appear at all. Such a simplifi- 
cation originates from a simpler scalar vertex in the case 
of a polarization operator. 

Nevertheless, the three methods discussed in the present 
Letter result in different, and hence, unphysical values for 
the interaction correction unless the large-momentum cut- 
off is imposed in the form of Eq. (29), or similar. In that 
case all methods yield the same value given by Eq. (3). 

To summarize, the calculations presented above indi- 
cate that the effects of electron-electron interactions lead 
to finite though numerically very small corrections to the 
minimal conductivity. Finally, the calculations of the 
present Letter are performed in the limit of zero temper- 
ature and their validity implies that Uu) » ksT. 

* * * 
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Appendix 

After publication of our paper in Europhys. Lett. 83, 
17005 (2008), a preprint by IF. Hcrbut, V. Juricic, 
O. Vafck, and M.J. Case, "Comment on "Minimal 
conductivity in graphene: Interaction corrections 
and ultraviolet anomaly" by Mishchenko E. G.", 
appeared in arXiv:0809.0725. It was previously under 
consideration for publication in the Europhysics Letters. 
Below we present our reply: 

The Comment argues against the procedure imple- 
mented above, which leads to Eq. (3), and advocates di- 
mensional regularization scheme in support of the value, 
Eq. (2), obtained in Ref. [18]. Yet, the Comment fails to 
offer a consistent resolution of the issue. Indeed, follow- 
ing our suggestion to utilize the charge conservation law, 
Eq. (5), the authors of the Comment analyzed the deriva- 
tion of the interaction corrections from the polarization 
operator and reported that the dimensional regularization 
yielded, C = (11 — 3tt)/6 ~ 0.26, the value different from 
their Eq. (2). (Note that this value coincides with akin 
given by Eq. (27) of the present paper when I is calcu- 
lated with the help of the hard cut-off.) Addressing this 
discrepancy, the authors of the Comment conclude only 
that, "The origin of this non- uniqueness is unclear at the 
moment, but we suspect that it may be the non-gauge in- 
variant contribution to the conductivity which is not prop- 
erly treated within the density polarization approach." 

Citing some unidentified contribution does not add clar- 
ity to the discussion. The equivalence between the Kubo 
and the density polarization approaches in the calculation 
of the homogeneous conductivity is ensured by the charge 
conservation law (the Ward identity) . It is surprising that 
the authors end the discussion with the above statement 
and do not even attempt to find out what happens to the 
charge conservation in their calculations. It is thus diffi- 
cult to conclude that the authors of the Comment were 
able " to clarify and correct some of the statements made" 
in our paper. If anything, the credibility of their result, 
Eq. (2), is even more questionable as it is now clear that 
this result is based on the scheme that yields values which 
vary depending on the method used. 

Interestingly, recent measurements of dynamic conduc- 
tivity [22] show a — <7o(1.01±0.04) over visible frequencies 
range and thus point towards smaller values of interaction 
corrections. 
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